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The scaling behavior of the current-voltage characteristics of chiral and gauge glass models of 
disordered superconductors, are studied numerically, in two dimensions. For both models , the linear 
resistance is nonzero at finite temperatures and the scaling analysis of the nonlinear resistivity is 
consistent with a phase transition at T = temperature characterized by a diverging correlation 
length £ oc t~ Vt and thermal critical exponent vt- The values of vt , however, are found to be 
different for the chiral and gauge glass models, suggesting different universality classes, in contrast 
to the result obtained recently in three dimensions. 
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Considerable attention has been paid recently to the 
nature of the glass phase and the universality class of the 
possible glass transitions in ceramic or granular high-T c 
superconductors. In the presence of an applied magnetic 
field, a vortex-glass phase transition with vanishing linear 
resistance in the low temnerature phase has been pre- 
dicted in three dimensionaxH while, in two dimensions, 
the linear resistance is nonzero at finite temperatures TL 
but there is, nevertheless, a zero-temperature transitional 
with a divergent correlation length £ cx T~ UT and critical 
exponent vt ~ 2 which determines the behavior of the 
current-voltage characteristics. The absence of a finite- 
temperature vortex-glass transition and current-voltage 
scaling has been .yenned in some two-dimensional super- 
conducting filmsu with a critical exponent vt in agree- 
ment with theory. In the theoretical studies, a gauge 
glass model has been widely used which is believed to 
be in the same universality class as the vortex glass. In 
zero external field, however, only the standard supercon- 
ducting transition takes place. In contrast, in d-wave 
superconductors materials containing " ir " junctionsu a 
chiral-glass phase has been predicted even at zero exter- 
nal field in three dimensions based on numerical studies 
of the. chiral glass modeloQ, or alternatively, and XY spin 
glassEl. The chiral-glass order parameter is the chirality, 
which represents the direction of local current loops in- 
troduced by frustration effects. As for the gauge glass, in 
two dimensions, tie .chiral-glass transition only occurs at 
zero temperature&oEj. It is well known that the chiral 
glass model has an additional reflection symmetry which 
gives rise to quenched in vortices in the proposed chi- 
ral glass phaseo. In spite of this, a recent study of the 
current-voltage characteristics of the chical and gauge 
glass models in the vortex representationEil found that 
they exhibit, within the numerical accuracy, the same 
critical exponents in three dimensions, suggesting a com- 
mon universality class. This would imply that resistiv- 
ity measurements alone are not able to distinguish be- 



tween chiral and vortex glass states. However, in two 
dimensions the problem has not yet been analyzed in de- 
tail although several studies have already shown that the 
zero-temperature chiral glass transition is characterized 
by two different divergent correlation lengthsHtrtd, £ c and 
£ s , with different critical exponents, v c ~ 2 and v s ~ 1 , 
describing chiral glass (or vortex glass order) and phase 
glass order, respectively. In particular the value of v c 
for the chiral glass turns put to agree with the exponent 
for the gauge glass models. The question then arises as 
to what correlation length is actually probed by nonlin- 
ear resistance measurements since this will determine the 
current-voltage scaling and the resulting behavior could 
be either consistent with vt — v Cl the same as the gauge 
glass model, or vt — v s which could serve to identify the 
zero-temperature chiral glass transition. 

In this work we present a numerical study of the 
current-voltage characteristics of chiral and gauge glass 
models, in two dimensions, in a representation in terms 
of the phases of the local superconducting order parame- 
ter. For both models, we find that the linear resistance is 
nonzero at finite temperatures and a scaling analysis of 
the nonlinear resistance is consistent with a phase tran- 
sition at T = temperature characterized by a diverging 
correlation length £ oc T~ UT and thermal critical expo- 
nent vt, in agreement with previous work. The values 
of vt , however, are found to be different for the chiral 
and gauge glass models, suggesting that measurements 
of nonlinear resistance probe mainly the phase correla- 
tion length and the models are in different universality 
classes, in contrast to the result obtained in three di- 
mensions in the vortex representation^. Thus, in two 
dimensions, measurements of nonlinear resistance could, 
in principle, be used to identify a possible chiral glass in 
two-dimensional ceramic superconductors. 

The chiral glass and gauge glass models can be de- 
scribed by the same Hamiltonian 
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where 9i is the phase of the superconducting order pa- 
rameter of a "grain "at site i of a regular lattice, J > is 
a constant Josephson coupling and screeruas-effects have 
been ignored . In the gauge-glass modelHEil, Aij rep- 
resents a quenched line-integral of the vector potential 
which is taken to be uniformly distributed in the inter- 
val [0, 2tt) , representing the combined effect of disorder 
and the external magnetic field , while in the chiral-glass 
modeO has a binary distribution, or ir , with equal 
probability, which may represent the phase shift across 
Josephson junctions in models of d-wave ceramic super- 
conductors even the absence of magnetic fieldliJ. Alter- 
natively, the chiral glass n^df-i^ just another represen- 
tation of the XY spin glassOu'Lrliill with random couplings 

J ij — i J . 

To study the current-voltage characteristics of disor- 
dered superconductors described by the Hamiltonian of 
Eq. (|l|), we assume a resistively shunted Josephson- 
junctipD (RSJ) model for the current flow between. 
grainsLLj and use an overdamped Langevin dynamicsEEl 
to simulate the noncquilibrium behavior. The Langevin 
equations can be written as 
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where I ext is the external current, rjij represents Gaus- 
sian thermal fluctuations satisfying 



< mj(t) >= 



(3) 



< %(*) Vkl(t') >= 2 -^-8 ljM 8{t ~ (4) 

and a capacitance to the ground C Q is allowed, in addi- 
tion to the shunt resistance R Q , in order to facilitate the 
numerical integratiorJij. We use units where fi/2e = 1 , 
R a = 1 , J = 1 and set the parameter JR 2 C = 0.5 
in the simulations, corresponding to the overdamped 
regime. The above equations were integrated numeri- 
cally using, typically, a time step 8t = 0.02 — 0.05 r 
(t = 1/R J ), time averages computed with 2 — 4 x 10 5 
time steps and the results averaged over 5 — 10 different 
realizations of the disorder. To determine the nonlinear 
resistivity (or resistance in two dimensions), p n i = E/J , 
an external current / is injected uniformly with density 
J = I / L along one edge of a square lattice of size L x L 
and extracted at the opposite one, with periodic bound- 
ary conditions in the transverse direction. The average 
voltage drop V across the system is computed as 
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and the average electric field by E = V/L . We have 
also computed the linear resistance, Rl = limj_>o E/J , 



without finite current effects, directly from the long-time 
equilibrium fluctuations of the phase difference across the 
system A6(t) = Eii( e U ~ °l,j)/L as 
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which can be obtained from Kubo formula of equi- 
librium voltage-voltage fluctuations, Rl — J dt < 
V(t)V(0) > , using the Josephson relation V = d6/dt 
. Lattices of sizes L — 16 , 24 and 34 were used in the 
simulations with the main results obtained for the largest 
system size. 
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FIG. 1. (a) Nonlinear resistance E/J as a function of tem- 
perature T for the gauge glass model. (b)Arrhenius plot for 
the temperature dependence of the linear resistance Rl- 

The nonlinear resistance E/J as a function of current 
density J and temperature T for the gauge glass is shown 
in Fig. |l|a and an Arrhenius plot of the linear resistance 
R L in Fig.|l|b. They are consistent with previous results 
obtained in smaller systemstl. The corresponding results 
for the chiral glass are shown Figs. |^a and |^b . The data 
for both models show the exgected behavior for a T = 
superconducting transitionoH. In Figs. [l]a and ^a, the 
ratio E/J tends to a finite value for small J, correspond- 
ing to the linear resistance Rl, which depends strongly 
on the temperature. This is confirmed in Figs. [IJd and 
||b where the linear resistance obtained at J = from 
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Eq. (JsJ) is consistent with an exponential decrease with 
temperature indicating a finite energy barrier for vortex 
motion. For increasing J , there is a smooth crossover 
to nonlinear behavior that appears at smaller currents 
for decreasing temperatures. If one assumes a T = 
transition with a power-law divergent correlation length 
£r oc T~ Vt and since the external current density J intro- 
duces and additional length scale I ~ kTj J, the behavior 
of the nonlinear resistivity normalized to Rl can be cast 
into the scaling fornxm 
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where g is a scaling function and g(0) = 1, which contains 
a single parameter, the critical thermal exponent vt ■ 
This scaling form implies that the characteristics current 
density J n i at which nonlinear behavior is expected to 
set in decreases with temperature as a power law J„i oc 
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from a fixed value, chosen to be 2 . The behavior is 
consistent with a power-law behavior and the slope of 
the loglog plot provides a direct estimate of vt = 2.2(2) 
. A scaling plot according to Eq. (Q) for the gauge glass 
is shown Fig. obtained by adjusting the parameter 
vt so that a best data collapse is obtained. The data 
collapse supports the scaling behavior of Eq. (|?|) and 
provides an independent estimate of vt = 2.4 . From the 
two independent estimates we finally get v T 9 = 2.3(2), a 
result consistent with a similar analysis of the nonlinear 
resistivity obtained frorn-,the coulomb-gas representation 
of the gauge glass modeB. Following the same procedure 
for the chiral glass model, we obtain vt = 0.9(2) from 
the loglog plot in Fig. ||a and vt — 1.2 for the best data 
collapse in Fig. ^b, giving a final estimate v T a — 1.1(2) . 
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FIG. 3. (a) Crossover current density J n ; as a function 
of temperature for the gauge glass model, (b) Scaling plot 
E/JR L x J/T 1+VT for v T = 2.35. 



FIG. 2. (a) Nonlinear resistance E/J as a function of tem- 
perature T for the chiral glass model. (b)Arrhenius plot for 
the temperature dependence of the linear resistance Rl- 

We now proceed to verify the scaling hypothesis and 
obtain a numerical estimate of the critical exponent vt 
for the gauge and chiral glass models. Fig. |^a shows 
the temperature dependence of J n i for the gauge glass, 
defined as the value of J where Ej JRl starts to deviate 



The values of i/j> for the gauge glass and chiral glass 
models obtained from the above analysis are quite dif- 
ferent from each other suggesting different universality 
classes for the T — transition. This is consistent with 
the additional reflection symmetry property of the chiral 
glassErO, where changing 0, — > — 6i leaves the Hamilto- 
nian of Eq. ([!]) unchanged, whereas for the gauge glass 
there is only a continuQus^ymmetry In fact, other stud- 
ies of the chiral glassfflll3 find that the T = transi- 
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tion is characterized by two different correlation lengths 
with different critical exponents, £ c oc T~ Vc describing 
chiral glass order and £ s oc T~ Us describing phase glass 
(or alternatively, XY spin glass) order with the estimates 
v c ~ 2.0 and v s ~ 1.0 . Since the relevant length scale 
for phase slippage, which leads to voltage fluctuations, 
at finite temperatures is given by £ s we expect that the 
nonlinear resistivity scaling can be described by Eq. (0) 
with vt — v s which agrees with our numerical result for 
. A more complicated scaling analysis is also possible 
involving two correlation lengths but since £ s << £ c at 
low temperatures, £ s should dominate the resistive be- 
havior and a scaling analysis with a single length scale is 
a reasonable approximation. 
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the phase correlation length which has different critical 
exponents for the chiral and gauge glass and could, in 
principle, be, used to identify, as has been done for the 
vortex glassu, a possible chiral glass in two-dimensional 
ceramic superconductors. 
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FIG. 4. (a) Crossover current density J n ; as a function 
of temperature for the chiral glass model, (b) Scaling plot 
E/JRl x J/T 1+VT for v T = 1.2. 

The distinct critical exponents found for the chiral and 
gauge glass models in two dimensions is in sharp contrast 
with the result obtained for the-same models in three di- 
mensions by Wengel and Youngjij where the critical expo- 
nents at the finite temperature transition agree within er- 
rors suggesting a common universality class. This implies 
that in three dimensions, resistivity measurement alone 
cannot distinguish between chiral and gauge glass states. 
However, our results suggest that, at least in two dimen- 
sions, measurements of nonlinear resistivity probe mainly 
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